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■ We report multiwavelength observations of the soft X-ray transient (SXT) 

^ ; XTE Jl 118+480, which we observed with UKIRT, HST, RXTE, EUVE and many 

other instruments and facilities. Adding radio (Ryle Telescope, VLA), sub-millimeter 
(JCMT) and X-ray [Chandra and SAX) data from the literature, we assembled the 
most complete spectral energy distribution (SED) of this source yet published. We 
followed the evolution of this source for 1 year, including 6 observations performed 
during the outburst, and one observation at the end of the outburst. Because of 
Jjj \ XTE J1118+480's unusually high galactic latitude, it suffers from very low extinc- 

tion, and its SED is near-complete, including EUV (Extreme Ultraviolet) observations. 
XTE J1118+480 exhibits an unusually low low-hard state (estimated inner radius 
350i? s ) and a strong non-thermal contribution in the radio to optical domain, which 
is likely to be due to synchrotron emission. We discuss the interstellar column density 
and show that it is low, between 0.80 and 1.30 x 10 20 cm -2 . We analyse the evolution of 
the SED during the outburst, including the contributions from the companion star, the 
accretion disc, the outflow, and relating irradiation and variability of the source in dif- 
ferent bands to the SED. We find no significant spectral variability during the outburst 
evolution, consistent with the presence of a steady outflow. Analysis of its outflow to 
accretion energy ratio suggests that the microquasar XTE J1118+480 is analogous to 
radio-quiet quasars. This, combined with the inverted spectrum from radio to optical, 
makes XTE J1118+480 very similar to other microquasar sources, e.g. GRS 1915+105 
and GX 339-4 in their low/hard state. We model the high-energy emission with a hot 
disc model, and discuss different accretion models for XTE J1118+480's broad band 
spectrum. 

Key words: stars: individual: XTE J1118+480, X-rays: binaries, ultraviolet: stars, 
optical: stars, infrared: stars 
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1 INTRODUCTION 

Soft X-ray Transients (SXTs) , also called X-ray novae, are a 
class of low mass X-ray binaries (LMXBs) . Among this class 
of sources, mo re than 70% are thought to contain black holes 
iCharlesil99ct) . The compact object accretes matter through 
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an accretion disc from a low-mass star via Roche lobe over- 
flow. The history of these sources is characterised by long 
periods of quiescence, typically lasting decades, punctuated 
by very dramatic outbursts. SXTs are usually discovered in 
X-rays, but outbursts are visible at every wavelength and in 
particular are often accompanied by radio activity. A typical 
outburst is characterized by soft X-ray emission dominated 
by thermal emission from the hot inner accretion disc, and 
optical/UV emission produced by reprocessing of X-rays. 
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Figure 1. Multi-wavelength lightcurve of XTE J1118+480 
in 2000: the 6 epochs of our observations are indicated by the 
HST points. Top panel shows RXTE/ASM 2-10 keV count 
rates in two-day averages. Periods with too few dwells for 
a precise measurement have been excluded for clarity. Sec- 
ond panel is HST/STIS far-UV fluxes derived from E140M 
spectroscopy and an upper limit from the late G140L ob- 
servation. Middle panel shows VSNET optical measurements 
( http://www.kusastro.kyoto-u.ac.jp/vsnet/Summary/jlll8.html I 
Open circles are two-day averages of visual estimates. Solid 
circles are individual electronic measurements. Triangles 
are pre-discovery photographic measurements. Arrows arc 
upper limits from all of these sources. The fourth panel 
shows K band fluxes measured by UKIRT. The lower 
panel shows 15 GHz radio data from the Ryle Telescope 
{http://www.mrao.cam.ac.uk/~guy/J1118-l-480/J1118480.list I. 
The last point is based on the average flux from 2000 August to 
November. 



One such source, XTE J1118+480, was discovered by 
RXTE on 2000 Marc h 29 at the galactic co ordinates (l,b) 
= (157.62°, +62. 32°) iRemillard et al.ll2000t) as a weak (39 
mCrab), slowly rising X-ray source, the post-analysis reveal- 
ing an outburst in January 2000, with a similar brightness. 
The outburst history of the source is shown in Fig. 

The optical counterpart of XTE J1118+480 in out- 
burst was a 13th magnitude star , coincident with a 18.8 
mag object in the USNO catalog llJemura et"afll200Cl) . Its 
optical spectrum w as typical of X-ray novae in outburst 
{Garcia et al.ll2000h . An infrared coun terpart was rapidly 
identified thereafter dChatv et al.ll2000h . Taking the X-ray 



and optical fluxes given above, it is clear that this sys- 
tem exhibits a very low X-ray to optical flux ratio o f 
5 (in /„, see iRemillard et alJ Eoool: lUemura et all l200(t l . 
compared to a typical value of 500 fo r SXTs in out- 
burst (see e.g. iTanaka fc Shibazakilll996t) . A weak 4.1 hr 
(0.17082 d) photom etric modulation was rapidly discovered 
jcook et al.l l200fj) . associated with the or bital period of 
the s ource, or possibly a superhump period dUemura et alJ 
2000). XTE J 11 18+480 exhibits therefore t he shortest or- 
bital p eriod among the black hole candidates. IWaener et alJ 
( 2001) took spectra of this source near quiescence and ob- 
served a period 0.5% smaller than during the outburst. 
They suggested that this difference was due to the pres- 
ence of super hump modulations during the outburst, later 
confirmed by IZurita et alJ l)2002f) . Flickering with an am- 
plitude of ~ 0.4 mag, and also a quasi-periodic oscillation 
(QP Q) at 10 s, wa s observed in the o ptical, in the UV 
ffiaswell et alJl200d) and also in X- ravs l| Revnivtsev et ail 
l2000h . withan evolving frequency dWood et al.l l200(t) . A 
faint radio counterpart was detected at 6.2 mjy, but no je t 
feature could be spatially resolved fsee lFender et ai]l200lT) . 
Optical observations of the source in quiescence led to the 
determination of a large value of the mass function, f(M) 
= 5.9 ± 0. 4 Mp), suggesting th at the compact object is a 
black hole JWaener et alJl200 it McClintock et al. 2001a). 

XTE J1118+480's location at an unusually high galac- 
tic latitude (b — +62°) in the direction to the Lockman 
Hole implies that there is a v ery low abso r ption along 
the line of sight of the source. iGarcia et al] J2000) esti- 
mated E(B-V) ~ 0.024, leading to a column density of 
Nh ~ 1.4 x 10 20 cm -2 . The column d ensity has been es - 
timated as N H ~ 0.75-1 x 10 2 cm" 2 jrlynes et al]l200d) . 
N H ~ 0.75-1.15 x 10 20 cm -2 jEsin et al]l200 ll) . Nh ~ 1- 
1.3 x 10 20 cm" 2 (McClintock et al. 2001b), N H ~ 0.7 - 
1.5 x 10 20 cm -2 (Frontera et al. 2001a), all nearly consistent 
with the mean value in t he direction of the Lockman Hole 
(0.5 - 1.5 x 10 20 c m" 2 . seelHvnes et al]l200fj) . From analysis 
of optical spectra. iDubus et al derived Nh ~ 1.77- 

4.47 x 10 20 cm" 2 . As we see, the exact column density is still 
a matter of debate, and will be one of the points discussed in 
this paper. McClintock et al. (2001a) estimated a distance of 
1. 8 ±0.6 kpc, con s istent with the value of ^ 1.5 kpc derived 
bv lUemura et al] feOOCh . 

We triggered our multi-epoch multi-wavelength override 
program with HST/ RXTE/VKIRT to get early observations 
of this system. We also requested Director's Discretionary 
EUVE observations because of the low value of interstel- 
lar ab sorption, and obtained the first EUV s pectrum of an 
SXT iMauche et al.ll200d . iHvnes et al.ll2000l) . We therefore 
have unprecedented broadband coverage, of more than 80% 
of the electromagnetic spectrum from the radio (A = 21 cm) 
to the 7-rays (180 keV) (see SED in Fig.0J. The analysis of 
one of the epochs (correspo nding in this paper to epoch 1, 
see Fig- El was described in lHvnes et al] i2000T) and we will 
briefly report here the main re sults (see also McClintock et 
al. 2001b and lEsin etafll200ll for an analysis of epoch 2). 
The corresponding SED suggested that the system was ex- 
hibiting a low-state mini-outburst, with the inner radius of 
the accretion disc estimated to be rather large, with a max- 
imum value at ~ 2000i? s (R s = 2S££. tue Schwarzschild 
radius for an object of mass M) demanded by the EUV and 
X-ray data. One of the most striking features was the strong 
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non-thermal (likely synchrotron) contribution in the optical 
and near-infrared (NIR) wavelengths, characteristic of out- 
flows seen in jet sources. Indeed, the SED shows a very flat 
spectrum from the UV to the NIR (~ 1000 - 50 000 A), 
suggesting that there is another source of NIR flux domi- 
nating any thermal disc emission, likely related to the radio 
emission, and therefore possibly of synchrotron origin. 

We followed the evolution of this system from the out- 
burst towards quiescence, to study the mechanisms under- 
lying the outburst. Preli minary results were published in 
IChatv et all feOOlaT) and IChatv et alJ fcOOlbl) . In this pa- 
per we will concentrate on the analysis and ev olution of the 
SED. A companion paper jHvnes et al.ll2003ft discusses the 
multiwavelength variability properties. We describe the mul- 
tiwavelength observations in Section |2] present the results 
in Section [21 and discuss them in Section [I] 

2 OBSERVATIONS 

Observational details for our VLA, Ryle telescope, JCMT, 
UKIRT, HST, EUVE, SAX, Chandra and RXTE data fol- 
low. 

2.1 Radio observations 

All the radio (VLA, Ryle Telescope) and sub-mm (JCMT) 
observations are taken from references reported in Tabled 

2.2 NIR UKIRT observations 

Near-infrared observations were carried out at the UKIRT 
3.8 m telescope using UFTI, IRC AM /TUF TI and also 
CGS4 ; some of which were already reported in lHvnes et alJ 
(2000). The log of the NIR observations is given in Table 
H The UFTI (UKIRT Fast-Track Imager) instrument is a 
cooled 1-2.5/^m camera with a 1024 x 1024 pixels 2 HgCdTe 
array. The plate scale is 0.091 arcsec per pixel, giving a field 
of view of 92 arcsec. The IRCAM/TUFTI instrument is a 
cooled 1-5/im camera with a 256 x 256 pixels 2 InSb array. 
The plate scale is 0.081 arcsec per pixel with a field of view 
of 20.8 arcseconds. 

Images were taken through the wide-band filters J98 
(A = 1.275/im; AA = 0.290/im), H98 (A = 1.670/xm; 
AA = 0.280/zm) and K98 (A = 2.205Aim; AA = 0.41/nm) 
with UFTI, and with all the above plus L'98 (A = 3.8/zm; 
AA = 0.6/Ltm) and M'98 (A = 4.675/zm; AA = 0.250mhi) 
with IRCAM/TUFTI. 

The exposure times range between 10 s and 60 s. The 
conditions were photometric. After taking each image of the 
object, an image of the sky was acquired, to allow subtrac- 
tion of the blank sky. The images were further treated by 
removal of the dark current, the flat field and the bright 
infrared sky. We also took a NIR K-band spectrum of this 
source using the CCS 4 instrument and a 0.6" slit on June, 
27.2 UT, which was featureless fsee lChatv et al.ll2003ft . 

2.3 Optical and Ultraviolet HST observations 

Hubble Space Telescope (HST) observations were performed 
with the Spac e Telescope Imaging Spectrograph (STIS; 
iLeithererlfeoOlll on the dates indicated in Table [3] These 



Epoch 


Date 


MJD 


Exp. time 


m 


8-9/04/00 


51642.603-51643.213 


19193.52 


13-14/04/00 


51647.821-51648.430 


19553.9 





16-19/04/00 


51651.407-51654.102 


80200.1 



Table 4. Log of the EUVE observations. 

The epoch, date, MJD, and exposure time in seconds are indi- 
cated. 

spanned the UV and optical bands at high and low resolution 
respectively using the E140M, E230M, G430L and G750L 
modes. For each visit, average calibrated spectra were ex- 
tracted from standard HST pipeline data products. There 
was useful coverage from 1150-10000 A, although the region 
from 1195-1260 A was co mpletely dominated by Lya ab- 
sorption and Nv emission (lHaswell et alj f2002l , so was ex- 
cluded from our spectral energy distributions. On the dates 
where more than one observation of the same wavelength 
range was taken, the spectra were averaged to increase the 
signal-to-noise ratio. The documented absolute calibration 
accuracy is 5% for the optical (CCD) modes, and 8% for the 
UV (MAMA) modes. The break between them is around 
3100A, with some overlap. This systematic uncertainty is 
larger than any statistical uncertainties. 

2.4 Extreme Ultraviolet EUVE observations 

Extreme Ultraviolet Explorer (EUVE) observations of 
XTE J1118+480 took place during 2000 April 8.10-8.71, 
13.32-13.93, and 1 6.91-19.60 UT. They are described in 
iHvnes et all fcOOdh The l og of the EUVE ob servations is 
reported in Tabled As in IHvnes et alJ (l2000tl . we consid- 
ered only the data between 70 and 120 A, since longer wave- 
lengths are heavily absorbed and dominated by noise. In the 
following plots showing data from EUVE, the vertical error 
bars are the 1 a errors from the photon statistics. 

2.5 X-ray SAX observations 

Beppo-SAX observed XTE J1118+480 4 times, on 2000 
April 14, May 4, June 26 and December 12. Here we show 
the observation taken on 2000 April 14-15, and reported in 
TableQ>] The details of this observation can be found in Fron- 
tera et al. (2001a). Since these data were initially corrected 
with an interstellar absorption of 1.5 x 10 20 cm -2 , we first 
uncorrected these data to get the observed flux, and then 
corrected them from the absorption with the desired value 
of Nh, as described in section |3~T1 

As shown in Figure [5] although there is a good agree- 
ment between SAX and Chandra data in the interval log v = 
16.6-17.1 and between SAX and _RJ\T"iJdata for the interval 
log v = 18.5-18.7, there is an inconsistency between the SAX 
and the EUVE observations for the interval log v — 16.5- 
16.7 (likely due to the model used in fitting the SAX data 
combined with the low response of the detector at those 
wavelengths) and in the interval \ogv — 17.0-18.5 between 
SAX and both RXTE and Chandra observations. This seems 
to be a calibration problem. 

Nonetheless, the SAX spectrum unambiguously shows 
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Epoch 


Date 


MJD 


Inst 


1.4 GHz 


8.3 GHz 


15 GHz 


23 GHz 


350 GHz 


Ref. 






30-31/03/00 


51634 


RT 






6.2±0.5 


_ 


_ 


Poolev fe Waldram 2000 


1 
1 




31/03/00 
02/04/00 
03/04/00 
05/04/00 


51635 
51637 
51638 
51640 


RT 
RT 
VLA 
RT 


■ 


■ 


7.8±0.35 
7.5±0.30 

8.7±0.3 


8.8 ± 0.3 




Dhawan et al. 2000 
Dhawan et al. 2000 
Dhawan et al. 2000 
Fender et al. 2001 




2-6 




16/03-25/06/00 


51620-720 


VLA 


2.6±0.4 


6.5±0.7 




9.3 ± 1.2 




Fender et al. 2001 




2-6 




16/03-25/06/00 


51620-720 


RT 






9.0±1.0 






Fender et al. 2001 








08-11/00 


51779-872 


RT 






0.15±0.17 






Pooley, priv. com. 
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30-31/05/00 


51695 


JCMT 










41±4 


Fender et al. 2001 








09/09/00 


51796 


JCMT 










< 21 


Fender et al. 2001 



Table 1. Log of radio observations. 

We tabulate here the epoch, date, Modified Julian Date, instrument and flux in mjy for every band. 



Epoch 


Date 


MJD 


Inst 


1 Z J 


H 


K 


L' M' 


m 


04/04/00 


51638.2 


UFTI 


12.12±0.02 


11.75±0.02 


11.06±0.02 










12-15/04/00 


51647.5 


Sternberg 


12.4±0.2 


11.9±0.1 


10.9±0.1 


9.2±0.1 


I 


2 ! 


1 


18/04/00 


51652.5 


TUFTI 


11.92±0.07 


11.43±0.06 


11.05±0.08 


9.71±0.14 9.38±0.42 








24/06/00 


51719.2 


UFTI 






11.512±0.004 






5 




26/06/00 


51721.2 


CGS4 






spectrum 






6 




15/07/00 


51740.2 


UFTI 






11.948±0.006 






7 




07/03/01 


51975 


UFTI 


17.41±0.05 17±0.05 16.72±0.05 


16.15±0.05 


15.77±0.05 





Table 2. Log of the infrared observations. 

The epoch, date, MJD, instrument and magnitudes for every filter are repo rted. The magnitudes of the fi rst epoch (corresponding to th e 
detection of the infrared counterpart) were reported in IChatv et al J 12000ft . and those from Sternberg in iTaranova fc Shenavrinl 12000 ) . 
The observations of epochs 5 and 7 are reported in more details in lGhatv et alJ 120031) . 



Epoch Date MJD Exp. time 

|~2~| 18/04/00 51652.5 27200 



Table 6. Log of the Chandra observations. 

The epoch, date, MJD and exposure time in seconds are indi- 
cated. For more details about the Chandra observations see Mc- 
Clintock et al. (2001b). 

2.6 X-ray Chandra observations 

There was only one observation by Chandra, taken on 2000, 
April 18, reported in Table|U Details are in McClintock et al. 
(2001b). Since these data were corrected with an interstellar 
absorption of 1.3 x 10 20 cm -2 , as for the SAX data, we first 
uncorrected these data to get the observed flux, and then 
corrected them for the absorption with the desired value of 
Nh, as described in paragraph 13. II 

2.7 X-ray RXTE observations and data analysis 

We observed XTE J1118+480 with the Rossi X-ray Timing 
Explorer (RXTE) Proportional Counter Array (PCA) and 
High-Energy Timing Experiment (HEXTE) at five epochs 
selected to coincide with the HST visits. The log of the 
R XTE observations is reported in Table|7| The method used 
for the analysis of data is the same than the one described in 




XTE/HEXTE 

XTE/PCA J| 
■i 



SAX 1 



Figure 2. Comparison of epoch 2 Spectral Energy Distribution 
as observed by the instruments EWE, Chandra, SAX and XTE. 
The SED is corrected with N H = 1.10 x 10 20 cm" 2 . 



the presence of a cut-off at an energy of log v ~ 19.5, as dis- 
cussed by Frontera et al. (2001a). Therefore in the following 
we will include the SAX data only at high energies (in the 
interval log v = 18.8-19.65) 



Multiwavelength observations of XTE J1118+480 5 



Epoch Date 


MJD 


E140M 


E230M 




G430L 


G750L 


Central wavelength (A): 


1425 


1978 


2707 


4300 


7751 


Wavelength rang 


;e (A): 


1 1 OQ 1 71 n 


1574-2382 2303-3111 


zyuu— 1 UU 


OZOO— lUZOO 




T 


08/04/00 


51642.7 


9150 


1300 


1200 


111 


180 




2 


18/04/00 


51652.6 


3000 


1000 


700 


1160 


150 




T 


28/04/00 


51663.3 


1620 


820 


700 


120 


150 




4 


28/05/00 


51692.8 


1800 


1000 


750 


120 


204 




T 


24/06/00 


51720.0 


1750 


1000 


700 


158 


216 








3000 












6 


08/07/00 


51733.9 


1700 


950 


750 


120 


221 




7 


28/09/00 


51815 


1250 


2050 


2050 


400 


407 



Table 3. Log of the optical and UV observations. 

The epoch, date, MJD, central wavelength and wavelength range are indicated, with in each case the exposure time in seconds. 



Epoch Date MJD 


LECS 


MECS 


HPGSPC 


PDS 


Energy bands (keV): 


0.12-4 


1.7-10 


7-29 


15-200 


|~2~| 14-15/04/00 51648-9 


21197 


29758 


42138 


20146 



Table 5. Log of the SAX observations. 

The epoch, date, MJD, energy bands and exposure time in seconds for the different instruments are indicated. For more details about 
the SAX observations see Frontera et al. (2001a). 



iHvnes et al.l i2000t) . therefore we will just give in the follow- 
ing the various parameters derived by analysing the different 
epochs. 

For the PCA, we used the "standard mode" data (128 
spectral channels, 16-second accumulations), selecting sub- 
intervals when the number of detectors on remained con- 
stant (about 90% of the total time). We similarly extracted 
256-channel spectral accumulations from the HEXTE sci- 
ence event (SE) data. A subset of PCA and HEXTE de- 
tector channels, corresponding typically to about 3-100 keV 
were used in our subsequent model fitting. Background rates 
for the PCA were estimated using the epoch-4 models, and 
response matrices were generated using the current calibra- 
tion files and response-matrix generation software, all from 
the "HEAsoft 5.1" release. 

The source intensity was typically in the 30-40 mCrab 
range for each epoch, with typical PCA count rates of 
~ 80— 120 cts/sec/PCU (source; the background is an addi- 
tional 30/s/PCU). The spectra thus derived were found to 
be hard, with photon power-law indi ces of about 1.8 ± 0.1. 
A th ermal Comptonization model iSunvaev fc Titarchukl 
1980) with t ~ 3 and T c ~ 30keV, also provided accept- 
able fits (in either the power-law or thermal Comptonization 
cases, typical x 2 P er degree of freedom of order unity were 
obtained) . There was no evidence (in terms of statistical im- 
provement to our fits) for a soft-excess component, thus we 
conclude, as have others, that the source remained in the 
"low/hard" spectral state throughout the outburst. In most 
cases, particularly 50133-01-03-00 for which the PCA expo- 
sure was about 10,000 seconds, there was a distinct positive 
residual corresponding to the 6.4 keV FeK resonance, thus 
a Gaussian line profile was included to refine the overall fit. 



Date 


MJD 


Exp. time 


Observation ID 


08/04/00 


51642.7 


3900 


50133-01-01-00 


18/04/00 


51652.6 


4700 


50133-01-02-00 


28/04/00 


51663.3 


10600 


50133-01-03-00 


28/05/00 


51692.8 


2400 


50133-01-04-00 


24/06/00 


51720.0 


5600 


50133-01-05-00 



Table 7. Log of the RXTE observations. 

We report here the epoch, date, MJD, exposure time in seconds 
for the PCA instrument and the observation ID. 



The energy coverage was ~ 3 to 26 keV with PCA and 11 
to 207 keV with HEXTE. We show in all the Figures only 
the data up to 80 keV (except for the first visit where we 
show them up to 120 keV), since at higher energies the noise 
dominates. 



2.8 Broadband SED 

For analysis, we separated the observations into the 7 epochs 
reported in Tables 00 Table summarizes the different 
facilities used at each epoch. For each epoch, observations 
in all bands are simultaneous or nearly-simultaneous. The 
broadband SEDs corresponding to all epochs are presented 
in Figures 0and0 In Figure0they are shown with different 
normalisations for more clarity. We overplot all epochs in 
Figure and enlargements in radio, NIR-UV and EUV- 
X-ray regions are shown respectively in Figures and 
In all these figures the data have been corrected with 
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Epoch Date 


MJD Radio 


UKIRT 


HST 


EUVE SAX 


Chandra XTE 


colour 




log(v): 


9.0-11.6 


13.78-14.48 


14.47-15.41 


16.38-16.61 16.61-19.68 


16.76-18.23 17.78-19.68 






1 


08/04/00 


51642.7 x 


X 


X 


X 


X 


dark blue 




T 


18/04/00 


51652.6 x 


X 


X 


X X 


X X 


light blue 




3 


28/04/00 


51663.3 




X 




X 


dark green 




4 


28/05/00 


51692.8 




X 




X 


light green 




5 


24/06/00 


51720.0 


X 


X 




X 


yellow 




6 


08/07/00 


51733.9 


X 


X 






orange 




7 


09/00-03/01 


51815-52467 


X 


X 






red 



Table 8. Log of the different epochs. 

The date and MJD are indicated. The "-" indicates that we used for the indicated epochs the observations of epoch 2, which is sensible 
since the SED did not change significantly in these wavelengths during 3 months. 



N h = 1.10 10' u cnT' 




N h = 1.1 lO^cm"' 



Figure 3. Spectral Energy Distributions of all the epochs, be- 
ginning from the top of the Figure. Colour-codings are as given in 
Tabled For easier reading, we multiplied epoch 1 by 10, and then 
divided epoch 2 by 10, epoch 3 by 100, epoch 4 by 10 3 , epoch 5 
and 7 by 10 4 and epoch 6 by 10 5 . 



N H = 1.1 x 10 20 cm" 2 corresponding to A v = 0.059, as will 



Figure 4. Spectral Energy Distribution of all the epochs over- 
plotted 



be described in section |3~T1 



3 RESULTS 

Before describing the results, we present in Section [3. li the 
method used to correct for interstellar absorption. 

3.1 Correction of interstellar absorption 

The interstellar absorption was corrected in the following 
way. First, we choose a value for the column density by 
fitting the UV and EUV fluxes and slopes, as will be dis- 
cussed in section 13.51 For example here we take Nh = 
1 x 10 20 cm" 2 . 

Using this value of column density, we then correct the 
EUVE data and SAX and Chandra X-ray dat a using the 
absorption cross sections of lRumph et al] (il994l) for HI, Hel 
and Hell with abundance ratios 1:0.1:0.01, typical of the 
diffuse interstellar medium. 




Figure 5. Radio Spectral Energy Distribution of all the epochs 
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N h = 1.1 lO^cm"' 




i 



Figure 6. NIR-UV Spectral Energy Distribution of all the 
epochs 



N„ = 11 10 cm ' 




Figure 7. EUV-X Spectral Energy Distribution of all the epochs 



Then, we assume that this inferred column density 
is interstellar, and we adopt an average gas-dust ratio of 



(N(H! + H 2 ) / 
iBohlin et all E 



N(H! + H 2 )/E(B - V)) = 5.8 



magnitude 



JBohlin et al.lll978ft . This leads to the value of E(B - V), 



equal to 0.017 in the example taken here. Taking the value 
of R v = A V /E (B — V) = 3.1 typi cal of the diffuse inter- 
stellar medium llCardelli et al.lll989l) . we get the value of A v 
(= 0.053 here). 

Finally, we correct the infrared (from UKIRT), optical 
and ultraviolet (HST) observa tions for inters t ellar absorp- 
tion with the extinction law o f lCardelli etakl dl989ll . using 
the derived value of absorption, leading to the inferred SED. 

We note that even in the FUV, reddening corrections 
are < 15 % (depending on E(B-V)) and that extinction curve 
and gas-to-dust scaling are not critical. Only in the EUV is 
the correction large. 



3.2 Geometrical parameters of the system 

We took a black hole mass of 7.2 ± 1.3 M (McClintock 
et al. 2001a), which corresponds to a Schwarzschild ra- 
dius of R s = 21 km. The mass ratio has been measured 
as Q = 1/q = = 27 ± 5 (an extreme value among 
SXTs, lOrosJ 12001)" therefore the mass of the donor star 
M 2 = 0.27 ± 0.05 Mrs,. The distance of the system has been 
determined as 1.71 ±0.05 kpc (see discussion in Section[!OJ 
and t he orbital period to P orb = 0.169937(1) days ~ 4.08 
hours jZurita et al.ll2002fl . 

This gives an o rbital separation a = 1.76 ± 0.1 x 10 9 
m (|Paczvnskilll97ll) . To derive the value of the outer ra- 
dius of the accretion disc, we take the intermediate value 
between the disc's tidal radius Rt and the 3:2 resonant ra- 
dius R23, since the source showed the presence of super- 
humps (see Section 0. The disc's tidal radius is taken as 
90% of the Roc he lobe radius, therefore R T /a = 0.58 ± 0.01 
"Eggletonlll983ft and the 3 :2 resonant radius is R 2 zfa ~ 0.47 
IWhitehurst fc Kindll99"l . Hence the outer radius we take 
is r ou t = 0.52a. The inclination of the system is chosen to 
be 70 ± 10° , consisten t with McClintock et al. (2001a, and 
also IZurita et al"ll2002ft . 

The inner radius of the accretion disc will be a free 
parameter, but will typically be between 300 and 450 R s - 
This corresponds to a low state, in contrast to the high state 
where the accretion disc extends very close to the compact 
object, ri n « 3 x R s , corresponding to the last stable orbit. 



3.3 Fit to the nearly-quiescent SED 

To characterize the nearly-quiescent system, we use only the 
epoch 7 data from UKIRT and HST. Since there is prob- 
ably still some contamination from the accretion disc, we 
take the combination of two emission models, one repre- 
senting the companion star and one representing the ac- 
cretion disc. The spectral type of the companion star has 
been determined a s K5-M1 V (McClintock et al. 2001a, 
IWagner et all 1200 if) . In our analysis, we therefore take a n 
Ml V star photosphere (T eff = 3400K. iBessell et aljll993l . 
The radius of the mass donor sta r can be estimated with 
R 2 /a = 0.15±0.01 jEggletonll983fl therefore R 2 = 0.39 R e . 
We take a black body model for the nearly-quiescent accre- 
tion disc. We therefore fit the epoch 7 UKIRT and HST 
data with three free parameters: the fractional contribution 
of the secondary to the total emission, the temperature of 
the accretion disc, and the distance of the system. We found 
that the best-fit parameters were a fractional contribution 
of 25 ± 2%, a remnant accretion disc at 6000 ± 50K and a 
distance of 1.71 ± 0.05 kpc with a reduced \ 2 = 1.3, and a 
tolerance of 10 -5 (see Figure|SJl. The fractional contribution 
of the secondary to the total emission which we derived is 
in agreement with th e contribution of ~ 28 — 36% found by 
IWagner et ail fcOOlj) and ~ 34 ± 8% by McClintock et al. 
(2001a) from optical spectroscopy. We also tried to replace 
the black body model by a power-law, but the fit was worse, 
showing that the synchrotron contribution had decreased 
between epochs 1-6 and 7. We finally varied the absorption 
through the column density interval allowed by the broad 
band SED (see section 1331 but this did not change substan- 
tially the best fit parameters. Therefore our simple model 
gives a result consistent with the other methods, and in the 
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N h = 1.1 lO^cm" 




Figure 8. Fit of the epoch 7 (the nearly-quiescent system), with 
a fractional contribution secondary /total emission of ~ 25%, a 
remnant accretion disc at ~ 6000-fC and a distance of ~ 1.7 kpc, 
taking into account the emission of a M1V star at T e ff = 3400K 
with a radius R2 = 0.39 Rq. 



following we will use the parameters derived from our anal- 
ysis, to add the contribution from the companion star to the 
overall fits at other epochs. 

It is interesting to point out that, by plotting the in- 
frared colour-magnitudes (observed on 2001, March, during 
the near-quiescence of the source, corresponding to epoch 
7) on a Hertzsprung-Russell diagram, we deri ved that the 
comp anion star had a spectral type of Ml V llChatv et alJ 
2003), in agreement with the spectroscopic results. On the 
other hand, the radius of the companion star that we used 
for the fit in Figure [8] is smaller than an isolated main se- 
quence star by a factor two. However, lHaswell et alJ ll2002l) 
showed that XTE J1118+480's companion star is an evolved 
star, consequently it may well have a different radius to an 
isolated main sequence star of the same type. 



3.4-1 Model of the accretion disk 

To fit the broad band emission due to the a ccretion disc 
we us e the simple parameterised model of iHvnes et alJ 
(2002). It is based on a combination of the classic viscously 



usly 

heated black body disc spectrum ( Shakura & Sunvacv 1973, 
iFrank et al J 11993) and th e modified temperature distribu- 
tion f or an irradiated disc (ICunninghamlll97f3 . IVrtilek et alJ 
1990). See these papers fo r derivations of the relevant tem- 
perature distributions, and lDubus et al.1 <ll999l) for a critique 
of the assumptions. 

The model spectrum is calculated by summing a se- 
ries of black bodies over radius. The local effective temper- 
ature of a disc annulus is determined by the emergent flux 
at that radius, such that T 4 ff oc fboi- The emergent flux 
is the sum of viscous energy release within that annulus, 
JVisc oc T 4 isc , and the X-rays reprocessed by the annulus, 
firr oc T^ t . Hence the effective temperature contains contri- 
butions from viscous heating (T oc R~ 3 ^) and irradiation 
(T oc R~ 3 / 7 ). The effective temperature profile of the disc 
is therefore represented by 

T e g(R) = Tvi sc (-R) + Tlt r (-R) . 

Both profiles are effectively controlled in the model by 
the temperature at the outer radius of the disc. The vis- 
cous temperature, T V i 3C (R out ) , is a free parameter, typically 
around 7500 K; the irradiation temperature, Ti rr (R ou t), we 
usually take to be Tin- (-Rout) = K (but see Section f3. 611 . 

The detailed fit is the sum of the secondary star (or 
nearly- quiescent system), the model of the accretion disc in 
the low state, and three different power-laws demanded by 
the data in the radio, NIR, and X-ray bands, with respective 
spectral indices 0.5, —0.15 and —0.8, taking the convention 



(1) 



As the figures show, these three power-laws are natu- 
ral fits to the SED. They have the respective expressions: 
/■radio = 7 x 10 ~ 31 x u °' 5 between 1 x 10 9 and 4 x 10 11 Hz; 
f UniT = 2.5 x 10~ 23 x v-°- 1B between 4 x 10 11 and 8 x 10 14 

-13 



Hz; and /„., 



X v 
1.2 x 10 



x v 



-0.8 



between 8 x 10 and 



10 1 



Hz. 



3.4 Epoch 2 

As we can see in Fig. [I] (which we will discuss in more de- 
tail in section l3~%ll . there were no gross changes in the SED 
between epochs 1 and 6, so in the following we will consider 
the epoch where we have the most data, i.e. epoch 2. Our 
analysis differs from previous ones (see Section because, 
in addition to the disc black body, we explicitely include con- 
tributions in different bands (for instance power-law spec- 
tra in radio, submm-NIR and X-rays), together with the 
emission in the NIR-optical domain of the nearly-quiescent 
system, using the results from Section 13.31 We are aware 
that this does not account for irradiation of the mass donor 
star, however the fits are just intended to be illustrative. 
XTE J1118+480 is a unique source in term of its low ab- 
sorption, allowing us to study its whole SED in great detail, 
derive some useful constraints, and analyse the evolution 
from the outburst to the near-quiescence. 



3.4-2 Examining the SED 

We show epoch 2 with the different models and fit in Figure 
[5] We see clearly two main characteristic s. Firstly the sourc e 
is exhibiting a very low low/hard state iHvnes et alJl20o"ol) : 
secondly there is a strong non-thermal contribution in the 
radio domain with an inverted spectrum, extending up to 
the UV wavelengths, and even in the X-rays. We will develop 
this later. 

In most of our fits the typical values of the inner disc 
temperature, Tin are between 20 and 30 eV and those of 
the inner disc radius R lu are between 300 and 4507? s (see 
Section f3. 51 below and Table [§Jl These values are consistent 
with those derived by McClintock et al . (2001b) but th e 
inner radii are higher than that derived bv lEsin et alJ l)200ll) , 
consistent with their model including a significant ADAF 
contribution to the EUV. 
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R in = 352 R s 
T viK = 7370 K 
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Figure 9. Spectral Energy Distribution of epoch 2 corrected 
with N H = 1.10 X 10 20 cm -2 . R in and T v i sc given in the Figure 
are those derived by our fits. The lower solid thick curve is the 
emission representing the near-quiescent system, the intermediate 
solid thin curve is the multicolour black body disc in the low state, 
and the upper dash dot thin curve is the multicolour black body 
disc in the high state. The straight dashed lines are the three 
different power-laws, as described in section [3.41 The solid curve 
following the data is the sum of the contributions from the near- 
quiescent system, the accretion disc in the low state, and the three 
power laws. 



3.5 The hydrogen column density, iVn 

We now try to better constrain the value of iVn by fitting the 
whole SED corrected with different values of A^h, with the 
inner radius of the accretion disc and the outer temperature 
as free parameters. All the other parameters are taken as 
described in previous sections. The results are reported in 
Table M The best fits had a reduced \ 2 of 54, with 728 
degrees of freedom. This high value of the reduced x 2 J us t 
shows how illustrative the fits are, because detailed spectral 
features such as the Balmer jump and the dip in the Chandra 
spectrum are not fitted. From the \ 2 values obtained and 
inspection of corresponding figures, reasonable fits appear 
to be with Ah between 0.8 and 1.3 x 10 20 cm" 2 . 

With Nh = 1.45 x 10 20 cm -2 (i? in = 313i? s ) no ac- 
cretion disc models are simultaneously consistent with the 
very different slopes in the UV and the EUV. With Ah = 
0.75 x 10 20 {Rin = 970i? s ) the presence of the X-ray power 
law makes it impossible to reconcile both high UV and low 
EUV fluxes. Ah = 1.1 x 10 20 cm" 2 is the only one of the 
values considered here which is consistent with both accre- 
tion disc models and the UV/EUV fluxes and slopes. We 
present the epoch 2 SED corrected with the column density 
AT H = 1.1 x 10 20 cm" 2 in Figure^ 

Therefore, from the results of the fits, and also by plot- 
ting these results on the observations, A^h = 0.8 — 1.3 x 
10 20 cm" 2 seems to be the preferred value. In the following 
we will consider the value of Ah ~ 1.1 X 10 20 cm" 2 as the 
most likely, and we will draw all the figures with this value. 



N H (xlO 20 ) 


Rin (Ra) 


T v i sc (K) 


T m (eV) 


x 2 


0.70* 


1369 


7373 


8.53 


65 


0.75* 


970 


7320 


10.97 


58 


0.80 


450 


7246 


19.32 


5 i 


0.85 


427 


7266 


20.15 


54 


0.90 


407 


7286 


20.94 


54 


0.95 


389 


7306 


21.72 


55 


1.00 


374 


7327 


22.44 


55 


1.05 


362 


7348 


23.06 


56 


1.10 


352 


7370 


23.62 


56 


1.15 


354 


7394 


23.60 


57 


1.20 


354 


7418 


23.67 


57 


1.25 


356 


7442 


23.65 


58 


1.30 


326 


7460 


25.27 


58 


1.35 


310 


7480 


26.37 


59 


1.40 


316 


7505 


26.08 


60 


1.45 


313 


7528 


26.35 


60 


1.50 


310 


7551 


26.62 


61 


1.55 


307 


7575 


26.90 


62 


1.60 


308 


7599 


26.92 


62 


1.65 


306 


7622 


27.13 


63 


1.75 


302 


7669 


27.57 


65 


1.90 


296 


7739 


28.25 


67 


2.05 


292 


7810 


28.80 


70 



Table 9. Results of different fits according to jVu- The iVn is 
fixed, i?i n and T v ; sc are free parameters. Although it is not an 
independent parameter in the fits, we give Ti n for information. 
The tolerance used in the fits was 1 X 10" 2 (except where an 
asterisk is written, where the tolerance was between 0.1 and 1). 
Reasonable fits appear to be with TVh between 0.8 and 1.3 X 
10 20 cm -2 , and by inspecting the figures 1.1 X 10 20 cm -2 seems 
to be the best value, and we will keep it for the rest of the paper 
(see text for a discussion of this). 



3.6 Irradiation 

In Figure 1101 we present the epoch 2 SED corrected with 
the column densities 1.10 x 10 20 cm" 2 , taking into account 
an irradiation with Tirr(-Rout) = 7370 K, as described in 
section 13.4.11 Figure [5] shows the same SED without any 
irradiation. We can see that the slopes in the optical/UV 
parts of the spectrum are not consistent with the presence 
of irradiation. However, the irradiation component in our fit 
presents lots of assumptions, mainly about the geometry of 
the irradiating/irradiated region. Therefore, in view of the 
uncertainties in the irradiation function, although it is not 
required by our characterization of the data (which is also 
very crude), some irradiation of the disc cannot be ruled 
out. 



3.7 Contributions to LboI 

The luminosity due to thermal-viscous dissipation in the 
accretion disc is given by integrating the luminosity aT 4 
between Ri n and Rout, for the two faces of the disc: 

L disc = 2 x 27rcrT 4 ut x R 3 out x (-L - - J— ) 

Adopting values typical of the fits in previous sections, 
i.e. R in = 352R S (corresponding to Tj n = 23.6 eV) and 
T vlsc = 7370 K, we obtain L di3C = 2.2 x 10 36 ergs" 1 . 
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N h = 1.1 10' 
R in = 352 R s 
T vise = 7370 K 
T, r , = 7370 K 




Figure 10. Spectral Energy Distribution of epoch 2 corrected 
with N H = 1.1 x 10 20 cm -2 , R in = 3b2R s and with an irradiation 



of T iTI (R out ) = 


7370 K (see 


section[3,4.1l. 






Radio 


NIR 


X-rays 


Frequency 


[9-11.6] 


[11.6-14.9] 


[14.9-19.3] 


Luminosity 


4.6 X 10 31 


5.3 x 10 34 


1.5 X 10 36 



Table 10. Contributions LboI m different parts of the SED 
(called radio, NIR and X-rays following the convention given in 
Section 13.4.11 . The frequency is given in log 10 (v/ Hz) and the 
luminosity in ergs -1 . 



In the three power laws corresponding respectively to 
the radio, NIR and X-ray domains, we integrate F v oc v a in 
the frequency range given in Table I1UI and multiply it by 
4-kD 2 (i.e. assuming isotropic emission). The correspond- 
ing luminosities are given in Tabic 1101 Therefore, assum- 
ing that Lx-rays is representative of the contribution from 
the corona, we have Ldisc > L COIona _, and we will discuss 
this result later in 14.3.11 This is consistent with the fact 
that in the case where all the gravitational power is dissi- 
pated in a static corona, nearly half of the coronal lumi- 
nosity intercepts the disc and is reprocessed/reflected so 
idisc ~ ^corona/2. Similarly if a fraction of the gravita- 
tional power is re leased in the disc (instead of the corona) 

corona. 

/2 jHaardt fc Maraschilll993Tl . 



3.8 The other epochs: evolution of the SED 

There were no gross changes in the SED between epochs 1 
and 6, so we overplot the fit to epoch 2 with the data from 
all epochs in Figure fTTI The SED did evolve a little during 
the outburst, and the best way to characterize this evolution 
is to quantify the change of the spectral index (defined as 
in Equation|TJ between optical (10 14 6 Hz) and X-ray (10 18 
Hz) domains. These domains are chosen because both show 
power-law spectra and we have simultaneous coverage dur- 
ing the six epochs of observations. This, reproduced in Fig- 
ure shows that the electron energy distribution remains 
the same during the whole outburst. From Figures |S| and 



[7| it appears that the fluxes decrease slightly, however the 
slopes do not change much. This suggests that the energy 
injected in the outflow decreases during the outburst. 

We tried to reproduce this evolution of the SED by mod- 
ifying some parameters of our simple model. One way to do 
this is by decreasing the outer (viscous) temperature. How- 
ever, changing this temperature modifies the slope of the UV 
part of the SED, and as we can see in Figure |S| this slope 
remains the same during the entire outburst. The other way 
to act on the multicolour black body disk would be to mod- 
ify the inner radius, but this will only modify the EUV part 
of the SED. Therefore, the only way to reproduce this SED 
evolution is to act on the power-laws, and we can indeed do 
this just by changing the constant of the power law (and 
not changing its exponent). This is consistent with a non- 
thermal contribution which would have a decreasing energy 
during the outburst, but the same particle energy distribu- 
tion. It is also suggestive of an outflow constrained to the 
central part of the accretion disc, whose size would decrease 
during the outburst, which could be correlated with the in- 
creas e in QPO frequency during the outburst iWood et alJ 

mi. 

This evolution is very interesting in the sense that it is 
different from what we see in other SXTs, and p articularly in 
other states. For instance, in XTE J1859+226. iHvnes et alJ 
(2002) observed a change in the NIR-optical SED that they 
modelled as a change in viscous and irradiation temperature 
(as described in section 13. 4.111 . In contrast, the evolution 
in XTE J1118+480 is mostly due to a change in the non- 
thermal (outflow) emission, compared to the thermal (disc) 
emission. 

The lack of pronounced evolution of the SED dur- 
ing 3 months is also reminiscent of the behaviour of jet 
sources such as GRS 1915+105 and GX 339-4, in the so- 
called "plateau state", where only small changes occur in 
the lightcurves. This, again, supports the idea of a steady 
outflow emanating from the source XTE J 11 18+480, as we 
will discuss in Section ^, 21 



4 DISCUSSION 

4.1 Variability and non-thermal contribution 

A comprehensive discussion about the short time-scale vari- 
ability of XTE J1118+480 during all the epochs of obser- 
vations in the different wavebands (in c luding NIR, UV and 
X-rays) can be found in lHvnes et alJ (|2003j). Here we just 
summarize the facts most relevant to our current analysis 
by combining our results on the SED and their results on 
the variability. IHvnes et alJ <|2003T) show that the XTE data 
exhibit a Poisson subtracted fractional RMS of about 35- 
37% for the first observations decreasing to 26-28% for the 
later ones. The UV variability is weaker than the optical or 
NIR, and stronger at longer wavelengths, with typical rms 
variability ~ 3 — 4 percent in the far-UV and ~ 4 — 5 percent 
in the near-UV. This suggests that at least some of the UV 
variability may be associated with the non-thermal compo- 
nent. 

The ~ 10 s QPO was seen in the optical, UV and X- 
rays, implying a common origin fo r this QPO through out the 
spectrum. The sampling time of IHvnes et al] (l2003lV s NIR 
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Figure 11. Spectral Energy Distribution: The fluxes are cor- 
rected with Nh = 1-1 x 10 20 cm -2 . The overall fit is a multi- 
colour black body disc model with an outer disc temperature of 
7370 K and inner disc radius of 352i?, s . Straight lines: different 
po wer laws, with spe ctral indices of 0.5, —0.15 and —0.8 (similar 
to lHvnes et al.ll2000ft . The lower curve corresponds to the black- 
body emission from the companion star. 
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Figure 12. Evolution during the outburst of the spectral index 
linking the optical (10 14,6 Hz) and X-ray (10 18 Hz) domains. 

observations was insufficient to search for the QPO, however 
they did detect flickering at NIR wavelengths, of larger am- 
plitude (~ 0.8 mag) than in the optical (~ 0.4 mag). This 
NIR variability, up to 50 %, is consisten t with a dominant 
non-thermal emission l|Hvnes et alj Eo03). Furthermore, the 
featureless NIR spectrum we took in 2000 June (see data 
plotted in yellow in Figure|3J is also consistent with the fact 
that the disc is not the only source of emission in this part 
of the spectrum (f or more details concer ning the NIR see 
IChatv et ai]l2001bl and lChatv et alJl2003l) . Hence both the 
SED and temporal behaviour suggest a strong non-thermal 
(likely synchrotron) emission at radio-UV wavelengths. Fi- 
nally, it also shows that the cut-off frequency characteristic 
of the synchrotron radiation might be in the optical domain 
or even at higher frequencies. 



4.2 How common is XTE J1118+480? 

4-2.1 SED: Comparison with other galactic sources 

Several aspects of the broad band spectrum of this source 
appear to be very similar to those of the other well stud- 
ied black hole binaries in the low state. The hard X-ray 
spectrum with a spectral index a ~ —0.8 and a cut- 
off around 100 keV is typical of low/hard state spectra 
observed in C ygnus X-l (e.g. Fron tera et al. 2001b) or 
GX 339-4 (e.g. lzdziarski et al.lll99l) . XTE J1118+480's X- 
ray PDS spe ctrum is very similar t o those of other black 
hole sources ijRevnivtsev et al.ll200fl() . However, the charac- 
teristic features are significantly shifted towards lower fre- 
quencies relative to other black hole sourc es with similar 
masses (cf Figure 2 oflRevniytsev et alJ200ol and Figure 1 of 
ISunvaev fc Revnivtsevll2^0ol) ~This suggests that the X-ray 
emitting region is larger in XTE J1118+480. Therefore, this 
hard X-ray variability, which is slower than in other SXTs, 
seems in turn to be consistent with a particularly large disc 
inner radius in this source. However, such "slow and ubiq- 
uitous" QPOs have already been seen in other black hole 
candidates in the low state. For instance, QPOs at P ~ 20 s 
have been observed both in X -ray and optical in the source 
GX 339-4 jMotch et alJll983l) . 

Another significant difference with other hard state 
sources is the low inner disc temperature inferred from 
our data. In the case of Cygnus X-l in the hard state 
the inner temperature is k Tj n ~ 150 eV (see e.g. 
Balucins ka-Church et al ] |l995l : Frontera et al. 2001b) con- 
trasting with our derived value kT m ~ 25 eV (Table |5J 
for XTE J1118 +480. This is cons istent with a "very low" 
low/hard state dHvnes et alj f2000). According to the mul- 
ticolour black body disc model, this lower temperature im- 
plies a larger disc truncation radius than in typical sources 
(a few hundreds R s instead of a few tens). This can also sim- 
ply be expressed in terms of ri„/r ou t: this ratio is 4-7 x 10~ 5 
for a typical SXT like XTE J1859+226 . which we also ob - 
served with similar intensive coverage llHvnes et al.ll2002T) . 
The ri„/r out ratio lies between 7 x 10 -3 and 2 x 10 -2 in the 
case of XTE J1118+480 (depending mainly on the value of 
the inner radius as derived with our fits). We also point 
out that the better energy coverage of the data available 
for XTE J1118+480 helped us getting accurate parameters, 
since the inner disc temperature of hard state sources is usu- 
ally poorly determined due to strong interstellar absorption. 

Finally, XTE J1118+480 differs from other transient 
sour ces by its exceptionally lo w X-ray to optical flux ra- 
tio jTanaka fc S hibazakJ Il996l) as well as its plateau-like 
light curve which contrasts with the exponential decay 
observed in many transients. However, once again in or- 
der to be consistent we have to compare XTE J1118+480 
with other sources in the low/hard state, and when this 
is done it becomes less peculiar. For instance, the X-ray 
to optical flu x ratio in the cas e of GX 339-4 can be as 
low as 2.5-3 jMotch et alJll98St . very similar to the value 
of 5 in the case of XTE J1118+480. Also, its plateau-like 
ligh tcurve is common in the je t source GRS 1915+105 (see 
e.g. Tanaka & Shibazaki 1996). 

Therefore, the main difference with other sources is 
the prominence of the non-thermal contribution. Indeed, 
the radio-NIR spectrum presents an inverted spectrum, 
typical of non-thermal optically thick synchrotron emis- 
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sion. This combined with a lack of pronounced evolu- 
tion of the SED strongly represents the signature of a 
jet. As a counterpoint, the radio spectral index during 
the 1999 outburst of XTE J1859+226 was negative most 
of the time, and the interval when it was positive was 
less than 4 days: its maximum value was 0.143 ± 0.180 
jBrocksopp et al]l2002h . c.f. the steady 0.5 value in the 2000 
outburst of XTE J1118+480. The radio-NIR component is 
much stronger in XTE J1118+480, exhibiting a low/hard 
state: in contrast the radio-NIR component was relatively 
weak in XTE J1859+226's high/soft state outburst. This is 
consi stent with the presence of an outflow (see e.g. iFendeil 
1200 lh . 

This signature of a jet has been observed many times in 
GRS 1915+105 ( see recent multiwa velength observation s by 
Ueda et al.ll2002l) Cygnus X-l (s ee IStirling et al.ll200ll and 
Fender et alJl2000h or GX 339-4 JCorbel et alJl200d) . It is 
usually correlated with the appearance of the low-hard X- 
ray state spectrum, which is believed to be the result of a 
coupling bet ween the Com pto nising corona and a compact 
jet (see e.g. iFendejEoOll and ICorbel et al]l2000h . We will 
discuss this in more detail in the Sections 14.2.21 and !4. 3. 21 

4-2.2 Comparison with galactic and extragalactic jet 
sources 

The strong non-thermal contribution attributed to an out- 
flow, implies that XTE J 11 18+480 is a microquasar. We 
therefore compare it to GRS 1915+105, the archetype of the 
microquasars. We can also compare XT E J1118+48 's SED 
with typical SEDs of qu asars given i n lElvis et alJ l)l994l) . 
following the analysis of lUeda et al.l <|2002T) in comparing 
GRS 1915+105 with typical quasars. 

To compare the ratio of energy of the outflow with the 
accretion energy, we need estimates of both. The outflow en- 
ergy can be estimated by taking the value vF v in the radio 
band (y = 10 10 Hz) for XTE J1118+480, GRS 1915+105 
and quasars. However, for the measure of the accretion 
power, we have to take the wavelength where an optically- 
thick thermal emission from the accretion disk presents a 
peak, therefore in the X-rays (v — 10 18 Hz) in the case 
of GRS 1915+105, in the EUV (y = 10 16 Hz) in the case 
of XTE J1118+480 and in the UV (u = 10 152 Hz) in the 
case of quasars. Then we compare thereafter the ratio of 
vF v between the radio and X-rays for XTE Jl 118+480 and 
GRS 1915+105 with the ratio between radio and UV for 
quasars. 

This outflow/accretion ratio is typically 10 -6 for 
XTE J1118+480, when it is between 10" 7 and 10~ 5 
for GRS 1915+105 respectively when this source is in a 
plateau state or exhibiting large radio flares. The value 
in XTE J1118+480 is intermediate between the two dif- 
ferent values exhibited by GRS 1915+105. The main dif- 
ference is that XTE J1118+480's ratio is stable, instead 
GRS 1915+105's ratio is varying on shorter timescales be- 
tween the two extreme values. This suggests that even if 
the contribution of the outflow to the total energy budget 
in the two sources is comparable, there is a large differ- 
ence in their behaviour, which is probably due to the fact 
that XTE J1118+480 is exhibiting a continuous outflow, 
and instead GRS 1915+105 shows energetic but sporadic 
ejections on short timescales. Calculations of energy budget 



ta king into account d ifferent ejection behaviours are given 
in lChatv et all fcOOll) by comparing GRS 1915+105 and fi. 

For the quasars the outflow/accretion ratio is typ- 
ically 10~ 2 and 1 0~ 6 respectively for radio-loud and 
radio-quiet quasars llUeda et aljl2002f) . This suggests that 
XTE J1118+480 (and the microquasars) fall into the regime 
of radio-quiet quasars. The inverted radio spectrum resem- 
bles that of radio-quiet quasars, and the analogue of the 
mm-break, around 10 13 Hz, exhibited by the radio-quiet 
quasars and related to self-absorption, seems to be the 
sub-mm excess, around 10 12 Hz, in XTE J1118+480. On 
the other hand, if the emission from the outflow also con- 
tribute substantially to the high-energy domain, then the 
outflow/accretion ratio becomes even bigger, up to the typ- 
ical values of radio-loud quasars. If this analogy is real, then 
the intrinsic luminosity of the jet could dominate the ener- 
getics of XTE J1118+480, as in the case of the quasars (see 
TabletTTjl. 

Therefore, although the high energy spectrum of 
XTE J 11 18+480 is very similar to what is observed in 
Seyfert galaxies, its phenomenology differs widely from that 
of strong extragalactic jet sources such as blazars. This is 
because blazars are dominated by relativistic beaming phe- 
nomena while the inclination in the case of XTE J1118+480 
precludes this. 

4.3 Accretion models 

The low X-ray luminosity as well as the relatively large in- 
ner disc radius inferred from the comparisons with the mul- 
ticolor blackbody disc model suggests that either the inner 
part of the accretion flow (i.e. at radii below R[ n ) is radia- 
tively inefficient (efficiency L^ Q \/{M x c 2 ) ~ 5 x 10 -4 ), or 
that a large fraction of material is ejected in an outflow, or 
both. 

The ADAF mod el advocated by llchimarul lll977l) and 
iNaravan &: Yi lll994h explains low accretion efficiency as 
a consequence of the inner part of the flow advected 
into the black hole . The ADAF mo del was applied to 
XTE J1118+480 by lEsin et all fcOOll) who found a good 
fit to the optical-hard X-ray spectrum. Another possibil- 
ity could be that the accretion power is advected into a 
jet or outflow, i.e. that a significant fraction of the energy 
could be used to power the jet. I ndeed ADAF solutions 
imply the flow is not bound (c.f. iBlandford fc Begelmanl 
1999). This outflow is indeed inferred from the flat radio 
spectrum and analogy with other sources in the hard state, 
where evidence for jets exists, as we discussed in Section 
14.21 The most striking example of such sources is Cygnus 
X-l that presents radio and X-ray spectra very similar to 
that of XTE J1118+480 and where a m illi-arcsec radio jet 
could be resolved Stirling et al]|200ll) . However, since in 
XTE J1118+480 the fraction of the luminosity that can be 
firmly attributed to the jet represents only a few percent of 
the total luminosity, (see Section f3. 71 . the jet would have to 
be radiatively inefficient. 

Actually the standard ID ADAF solution is likely 
to be strongly affected by the presence of the outflow 
and, more generally by the development of convective 
instabilities that only a ppear in 2D solutions (see e.g. 
[ Abramowicz et alJ l2000t INaravan et all |2000L IStone et al. 
^999[ iLlmrnen^hchev et al] l200oT Touefeiert^^Gru^rrio"^ 
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2000). As a consequence of the large theoretical uncer- 
tainties, we prefer to illustrate the possible presence of a 
hot optically thin geometrically thick component at the 
center of the cold outer disc using a more phenomenological 
approach based on the constraints given by the observations 
and energy balance conditions as described below. 

4-3.1 The hot disc model 

We will now model the optical-hard X-ray spectrum, by as- 
suming that the hot component comprises a sphere with ra- 
dius equal to the cold disc truncation radius. The geometry 
is very similar to tha t of the sphere+disc geometry detailed 
in lDove et all dl997l) . 

To limit the number of parameters, and to avoid the 
uncertainties about the strength of the magnetic field, we 
further assume that the dominant cooling mechanism in the 
hot plasma is Compton cooling by the soft photons produced 
in the external disc. Other possible cooling mechanisms 
(bremsstrahlung, thermal synchrotron) are assumed to be 
negligible. This is a good approximation for luminous ac- 
creting black hole sources (' Wardzinski fc Zdz iarski 2000(). In 
the case of XTE J1118+480 however, cyclotron/synchrotron 
and bremsstrahlung could be important cooling mechanisms 
(c.f. Frontera et al. 2001a), but are not necessarily. We will 
discuss this further below, and will show how our assump- 
tion turns out to be consistent with the data. 

We therefore used the non-linear Monte-Carlo code 
of lMalzac fc Jourdainl ((2000) to compute the Comptonised 
spectrum emitted by the hot phase (the inner accretion flow) 
radiatively coupled in energy balance with the outer cold 
standard disc. In practice, we assume a homogeneous den- 
sity, dissipation and heating rate inside the hot sphere. For 
a fixed Thomson optical depth, the hot plasma temperature 
is computed by balancing the heating and cooling. Possible 
temperature gradients due to the inavoidable inhomogene- 
ity of the Compton cooling are accounted for by dividing 
the hot sphere in 10 homogeneous zones with equal volumes 
where the energy balance is computed locally. We neglect 
the effects of irradiation on the disc temperature profile i.e. 
we assume the standard viscously heated multi-colour black 
body disc (T cx R~ 3 ^ 4 ). The reprocessed emission is ac- 
counted for by re-emitting the absorbed energy, at the point 
it impinged on the external disc, with a black body spectrum 
at the local disc temperature. 

The escaping X-ray spectrum is controlled mainly by 
the ratio of the volume averaged dissipation rate in the hot 
phase (electron heating) to the soft flux from the disc that 
enters the cold phase (controlling the cooling). This ratio 
defining the energy balance of the hot component depends 
both on the assumed geometry and the fraction of radiated 
power dissipated in the hot phase: 

/ = L C /(L V + L c ) (2) 

where L c is the power radiated in the hot phase, L c ~ Lx- 
L v is the power viscously dissipated in the outer disc. The 
cold disc emission arises from both internal viscous dissi- 
pation L v but also from reprocessing of the Comptonised 
hard X-rays irradiating the external disc. For our assumed 
geometry approximately 1/3 of the hard X-ray luminosity 
is intercepted by the disc. The observed disc luminosity is 
thus : 



L disc = [(1 - /) + f/3] (L v + L c ) (3) 

and in the case of XTE J 11 18+480, / can be determined 
observationally: 

/ ~ + ^) ~ 0.33 (4) 

The fraction of disc luminosity due to reprocessing is 
small, ~ 3-2/ = 0.14 and irradiation is unlikely to affect 
the outer disc structure. In particular, the local tempera- 
ture depends only weakly on the total flux (T cx -F 1//4 ), thus 
for a local irradiating flux of 20 % of that of the local vis- 
cous flux, the disc temperature increases by less than 5 %. 
We note that this is perfectly in agreement with our ob- 
servations suggesting that effects of disc irradiation on the 
temperature profile should be weak (if any, see Section 
In the framework of the sphere+disc model, this "weak irra- 
diation" hypothesis is corroborated by the relatively low ob- 
served ratio of Comptonised to thermal emission. Therefore 
our simple multi-colour black body and power law fits are 
consistent with our more sophisticated model assumptions. 
We also note that this geometry predicts an amplitude for 
the reflection component R ~ 0.3 which is slightly above the 
upper limits obtained by Frontera et al. (2001a) for their fits 
assuming low met allicity abundances in the disc (but see also 
Mille r et al.ll2002T ). However, even if such low Z abundances 
mi ght be related to a poss ible halo origin, as argued for 
bv lWaener et al.l (1200111 and lMirabel et al] ]200ll) . we point 
out that the metallicity should be higher than expected for 
a halo object, s ince the mass donor star is probably a CNO 
processed core (IHaswell et al.l l 20021 . 

The spectral shape depends on 5 parameters: the disc 
outer temperature T out and radius Rout, the disc inner ra- 
dius Rin, and the fraction of accretion energy dissipated 
in the hot sphere to that dissipated in the accretion disc 
/, the hot sphere optical Thomson depth r with respect 
to the sphere radius. The unique broad-band coverage in 
XTE J1118+480 constrains all the parameters relatively 
well. 

In the simulation shown in Figure fTTH the cold disc 
parameters (and especially the locally emitted black body 
spectrum at its surface corresponding to the standard tem- 
perature dependence) have been set to values close to that 
assumed or derived from the fits of section 13.41 namely 
R out = 12 x 10 3 R s (80 % of the Roche lobe radius), 
Ri„ =300 R 3 , T ou t = 8000 K. The Thomson optical depth 
along the sphere radius was fixed at r = 1 as indicated by 
the fit of Frontera et al. (2001a) for a spherical geometry, 
and / = 0.33 as discussed above. We did not include in 
the model any absorption as seen at log(^) ~ 17 wh ich was 
attributed to a warm absorber bv lEsin et al] ]200ll) . 

The most important parameters for the X-ray emission 
are / and r. / controls the Compton y parameter of the 
plasma at equilibrium y cx tT c and thus the X-ray photon 
index F is very sensitive to /. With the / value indepen- 
dently provided by the observed luminosity of the cold and 
hot components, the model produces V = 1.82 (2.-10 keV), 
in agreement with the results from spectral fits of the SAX 
data (Frontera et al. 2001a). Since all the parameters deter- 
mining the X-ray slope are independently constrained, such 
an agreement is very remarkable. At fixed /, the parameter 
r controls mainly the temperature of the hot phase. From 
energy balance we get a volume averaged equilibrium tern- 
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Figure 13. Modelling the optical to hard X-ray spectrum. The 
hot disc scheme is illustrated here, where we compare the non- 
linear Monte-Carlo simulation for a central spherical hot plasma 
which radiates through Comptonisation of the soft photons emit- 
ted by the outer cold disc both through internal viscous dissipa- 
tion and reprocessing of the Comptonised hard X-rays (the later 
being fa 20 % of the disc emission). See section f4.3.1l 

perature of the sphere of kT e — 109 keV in agreement with 
the value inferred from the fits of SAX data (Frontera et 
al. 2001a). As a consequence, the shape of the high energy 
cut-off is well reproduced. 

Note that th is model differs from the ADAF model of 
lEsin e t al. ( 200JJ) who considered cyclo-synchrotron radia- 
tion as the main source of soft photons for Comptonisation 
and neglected the soft-photons from the disc. We did pre- 
cisely the opposite. In our model, if the bremsstrahlung and 
cyclo-synchrotron radiation were to constitute a significant 
additional source of soft photons, this would increase the 
cooling rate in the hot plasma, and, as a result we would 
get a much softer spectrum than was observed. The obser- 
vations would then require a different geometry, reducing 
the soft photon input from the outer disc. For example an 
oblate geometry for the central hot region could satisfy this 
requirement. In the extreme limit where the hot phase is ge- 
ometrically thin alo ng the disc axis w e would be in the situ- 
ation considered bv lEsin et alJ i200ll) . dominated by cyclo- 
synchrotron radiation. Note that the correlated variability 
between the optical and X-rays suggests a causal connec- 
tion, and that the synchrotron radiation emanating from 
the same region th at produces the X-r ays can contribute 
in the optical band jMerloni et ahfe oOO ) , but not necessar- 
ily as a significant source of soft cooling photons. Therefore 
this processus is not inconsistent with our model where we 
neglected the synchrotron emission. 

4-3.2 jet corona models 

If the X-ray emission is produced in the central part of the 
accretion flow, the ~ 0.1 Hz QPO obs erved both in the X - 
ray and optical is difficult to explain jMerloni et al-l boOO). 
On the other hand it may be explained simply by a flux 
modulation at the Kepler frequency of the disc at Ri„. This 
suggests that a non-negligible fraction of the flux would be 



produced in a transition region where both the hot plasma 
and cold thin disc coexist. 

This transition region could consist of a hot corona 
atop the external disc. The generation of magnetic field 
in the disc through magnetorotational ins tability and the 
subsequent buoyancy of th e magnetic field llTout fc Pringlj 
Il99l iMiller fc Stone] l200d) is generally invoked as an effi- 
cient mechanism for the transfer of magnetic energy from 
the disc to a corona where it is dissipated through magnetic 
flares. This scenario for energy dissipation in the corona was 
built to explain a strong coronal emission in the innermost 
part of the accretion flow. At the large distances we infer for 
the disc truncation radius, it is expected to be less efficient. 

Therefore, if Ri n is as large as we infer, alternative dis- 
sipation proce sses such as transport and heating through 
Alfven waves llTagger fc Pellaf j|l999l) or viscous dissip ation 
in a fast accreting corona iRozariska fc Czernvl 120001) may 
appear preferable. Actually, if the inner flow was very in- 
efficient (or even nonexistent), the X-ray emission could 
be fully dominated by the coronal emission. In this con- 
text, mildly relativistic ejections of coro nal plasma are likely 
to be important in hard state sources jBeloborodovl ll99SL 
iMalzac et all 1200 ll) . In particular it may explain the rela- 
tively hard X-ray spectrum and absence of reflection ob- 
served in XTE J1118+480 together with being consistent 
with the presence of an outflow. From the theoretical side 
relativistic outflow models often require an accretion disc 
corona as the place where the flow is initially powered (see 
iTagger fc Pelladll999l. iMerloni fc Fabianll200ll) . The accre- 
tion disc corona picture would be thus physically compatible 
with the presence of a jet inferred from the flat radio spec- 
trum. 

In accretion disc corona models a substantial part of 
the disc emission is due to reprocessing of the hard X-rays 
impinging on the disc. As the optical, near-UV and far-UV 
variability indicates that reprocessing is weaker than the ex - 
tremely strong synchrotron component jHvnes et alJl2003l) . 
this constrains the corona to be active only in the inner- 
most part of the disc. The X-ray luminosity should then 
come from a region of the disc forming a thin ring with 
inner radius Ri n and outer radius R c . 

One can get an estimate of R c by assuming that at dis- 
tances lower than R c , all the accretion power goes to the 
corona while at larger distance it is dissipated viscously in 
the disc. Since the estimated X-ray luminosity is about half 
of that of the disc and about half of the coronal emission 
is intercepted by the disc we can conclude that about half 
of the disc luminosity could be due to the reprocessing. As- 
suming that reprocessed luminosity is half of that of the disc 
one gets r c ~ 2R[ n . 

We point out that our results are perfectly consis- 
tent with the Accretio n-Ejection Instability (AEI) model 
( Tagger & Pellafj 119991) which predicts that Alfven waves 
efficiently power the corona and the jet in a range of dis- 
tance comprised between Ri n an d the corotation radius o f 
the Rossby vortex at a few _Ri n llVarniere fc Tagged 1200 j) . 
Particularly, from the value of the inner radius that we de- 
rive with o ur fits, the AEI predi cts a QPO which is well in 
agreement llVarniere et al]l2002l) with the values observed 
in UV, optical and X-rays during the outburst at Fourier 
frequencies ~ 0.1 Hz. The prediction on the QPO frequency 
evolution would also be consistent with the observations by 
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IWood et alJ i200Ct) . The AEI context, by assuming equipar- 
tition between magnetic field and gas pressure, would ex- 
clude a Schwarzschild black hole, and favor a Kerr black 
hole rotating at a spin between 0.90 and 0.99 (P. Varniere, 
priv. com.). This QPO may be explained simply by a flux 
mo dulation due to the kep lerian rotation of the disc or, in 
the lTagger fc Pellatl Jl999l) framework, the presence of the 
spiral wave (P. Varniere, priv. com.). 

We also note that in addition to the Comptonised emis- 
sion from the hot disc and/or corona, a synchrotron com- 
ponent from the jet could also contribute in the X-rays. 
Such a non-thermal component was seemingly detected i n 
GRS 1915+105 in a low/hard state dVadawale et al.ll200ll) . 
In fact, the whole X-ray emission could even arise from pure 
synch rotron from the jet as demonstrated bv lMarkoff et all 
(2001). This model attributes all the X-ray luminosity of 
XTE J1118+480 to synchrotron emission from the jet. A po- 
tential problem for this model is that it predicts essentially 
no reflection component. This makes it difficult to trans- 
pose to other similar jet sources such as Cygnus X-l where 
a significant reflection component is clearly observed and 
correlated with ot her X-ray spectral ch aracteristics such as 
the photon index dGilfanov et al.l lT999 i *l. However the pure 
synchrotron model cannot be formally ruled out, and we 
could also have both mechanisms acting at the same time. 



5 CONCLUSIONS 

We reported multiwavelength observations of 
XTE J 11 18+480 during its outburst, assembling the 
most complete spectral energy distribution (SED) of this 
source yet published, including our observations with 
UKIRT, HST, RXTE, EUVE, and adding observations 
from the literature: Ryle Telescope, VLA, JCMT, Chandra 
and SAX. We followed the source for 6 months, and show 
its evolution during the outburst. The main results of our 
broadband multi-epoch coverage are: 

(1) The source XTE J1118+480 was in a very low state 
throughout the outburst (estimated inner radius at 350i? s ). 

(2) The column density is low, between 0.80 and 1.30 x 
10 20 cm" 2 . 

(3) The accretion disc seems to be heated mainly by vis- 
cosity throughout the outburst, without a strong irradiation 
contribution. 

(4) It exhibited an inverted spectrum from radio to 
at least optical wavelengths, characteristic of a strong non- 
thermal (likely synchrotron) contribution, usually attributed 
to an outflow. 

(5) By examining the near-quiescent system we found 
that the best-fit parameters were a fractional mass donor 
star contribution of 25 ± 2%, a remnant accretion disc at 
6000 ± 50-ft: and a distance of 1.71 ± 0.05 kpc. 

(6) The quasi-absence of variation of the SED during 
3 months is consistent with a steady outflow emanating 
from the source, similarly to other "microquasars" in the 
"plateau" state. 

(7) The ratio between radio and EUV energy, indicative 
of the outflow to accretion energy ratio, suggests that the 
source falls into the regime of the radio-quiet quasars, and 
is also consistent with a steady outflow. 

(8) The small change of the SED in the optical-UV part, 



along with the constant power-law slopes during the whole 
outburst, is consistent with a continuous, gradual decrease 
of the outflow energy. 

(9) We modelled the emission from the optical to the 
hard X-rays with a hot disc model, showing that the high- 
energy part of the spectrum can emanate from the accretion 
flow. We therefore have to take into account the possibility 
that the high-energy emission from this source comes from 
i) Comptonization ii) pure synchrotron or iii) a mixture of 
both. 

Although this object exhibits peculiar characteristics, 
e.g. a large inner disc radius and a likely origin in the halo, 
some of its characteristics are very similar to other sources, 
particularly in the low-hard state and exhibiting jets or out- 
flows. This object is particularly important because its very 
low absorption reveals phenomena that until now were dif- 
ficult to study. XTE J1118+480 facilitates testing and re- 
finement of models for variability and emission in black hole 
accretors. 
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